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Abstract 

The structure and algorithms of the Monte-Carlo generator ELRADGEN 2.0 de- 
signed to simulate radiative events in polarized ep-scattering are presented. The 
full set of analytical expressions for the QED radiative corrections is presented and 
discussed in detail. Algorithmic improvements implemented to provide faster simu- 
lation of hard real photon events are described. Numerical tests show high quality 
of generation of photonic variables and radiatively corrected cross section. The com- 
parison of the elastic radiative tail simulated within the kinematical conditions of 
the BLAST experiment at MIT BATES shows a good agreement with experimental 
data. 
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1 Introduction 

The exclusive photon production in lepton-nucleon scattering is the routine 
experimental tool in investigating the hadronic structure. Depending on the 
design of experiments, the measurements of this process can give an access 
to the generalized parton distributions [T|E] or the generalized polarizabilities 
PP] . In some cases the exclusive photon production appears as a background 
effect to inelastic [5|6] or elastic [7] lepton nucleon scattering. The last scenario, 
i.e. the situation when the events with the real photon emission accompany 
the elastic electron-proton scattering is the most advanced due to the infrared 
problem, therefore it will be in our main focus. 

The set of processes contributed to the observed cross section in the next order 
of perturbation theory is referred to as the lowest order radiative corrections 
(RC). The basic contribution to the lowest order RC appears from the square 
of amplitude that only includes real photon emission from the lepton leg. This 
contribution contains the so-called large logarithm {i.e., the logarithm of the 
lepton mass) and normally is only held in the lowest order RC. 

In practice of data analysis, RC are calculated theoretically or their contribu- 
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tion to the observed cross sections (or asymmetries) are minimized by exper- 
imental methods. Due to finite detector resolution, a complete removal of the 
events with radiated hard photon(s) by pure experimental methods is not pos- 
sible. Furthermore, the contributions of additional virtual particles and soft 
photon emission cannot be removed in principle. The theoretical calculation 
provides with analytical expressions included the contributions of loops and 
photon emission which are infrared free after the procedure of the cancellation 
of the infrared divergence. The contribution of the hard photon radiation is 
presented in the form of integrals over photon phase space. Partly the inte- 
gration is performed analytically without additional simplifying assumptions 
or assumptions on specific functional forms describing hadronic structure. 

The pioneering approach for RC calculation in inelastic processes was sug- 
gested by Mo and Tsai in their seminal paper in 1969 [8] . They also developed 
the peaking approximation allowing for analytical estimating integrals over 
photon angles. The approximation is used in many data analysis, e.g. , in ref. 
[9] the electromagnetic RC in elastic ep-scattering was calculated in peaking 
approximation with taking into account the one-photon emission both from 
lepton and hadron legs. 

The Mo and Tsai approach requires involvement of the artificial parameter A 
separating the integration region over photon energy on parts with soft and 
hard photon contributions. To cancel infrared divergence analytically only 
leading terms in the expansion of the soft photon contribution over recipro- 
cal of the photon energy are kept. As a result, the final expressions contain 
undesired dependence on the artificial parameter. 

Bardin and Shumeiko developed the approach [TOj for exact separation and 
cancellation of infrared divergence when the final expressions for RC were com- 
pletely free from any artificial parameters like A. Using this approach, RC in 
polarized elastic ep-scattering withing QED theory has been calculated in refs. 
pT]|T2] . Basing on these calculations the FORTRAN code MASCARAD has 
been developed and successfully used for data processing of the relevant par- 
ity conservation experiments p^fTi] . Other approaches were also used for RC 
calculation in elastic ep-scattering. Thus, the total lowest order RC (both to 
lepton and hadron legs) was also calculated in [7] with soft photon approxi- 
mation and the method of electron structure functions suggested in the work 
[15j was also applied for estimation of RC to elastic ep scattering [T6|IT7] . 

The use of realistic detector geometry requires essentially complicated inte- 
gration over the real photon phase space. As a result, the researchers come to 
the necessity of using the Monte Carlo technique which constitutes a comple- 
mentary approach to the theoretical calculations of RC using respective codes 
such as MASCARAD. The Monte Carlo generators for simulation of radiative 
events have been developed for many specific processes and intensively used in 
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Fig. 1. Feynman graphs contributing to radiatively corrected cross sections of elastic 
lepton-nucleus scattering: Born (a), additional virtual particles (b,c) and real photon 
emission (d,e) contributions. 

data analysis. Thus, the Monte Carlo generator RADGEN [18] for simulation 
of radiative events in inclusive deep inelastic scattering has been developed 
on the basis of the FORTRAN code POLRAD [19j. The Monte Carlo gener- 
ator MERADGEN [20] for simulation of radiative events in M0ller scattering 
appeared on the base of FORTRAN code MERA [21]. 

In this paper we present and describe in detail the latest version 2.0 of the 
Monte Carlo generator ELRADGEN. The prototype of the code [22] dealt 
with the simulation of real hard photon emission as a background effect in the 
unpolarized elastic electron-proton scattering. The present version 2.0 is ex- 
tended on the initial polarized particles: longitudinally polarized electron and 
arbitrary polarized proton. The theoretical background for the developments 
is presented in ref. [TT] . 

The paper is organized as follows. Section [2] describes the kinematics of the 
investigated process and the generation method. The different contributions 
to the lowest order RC and the multi-soft photon emission are presented and 
discussed in Section [31 The brief structure of the code and the input-output 
datafiles are described in Section |H Test runs, comparison with MASCARAD, 
and numerical comparison of the simulated and measured cross sections of the 
radiative tail from elastic peak in the BLAST experiment are presented in Sec- 
tionO Conclusions and final remarks are given in Section|6l The four-momenta 
reconstruction formulae, explicit expressions for the lepton and target polar- 
ization vectors, some lengthy formulae for RC, and test outputs are given in 
Appendices. 



2 Kinematics and Method of Generation 

The lowest order (Born) (FigJT] (a)) as well as the additional virtual particle 
(FigH] (b,c)) contributions to the polarized elastic lepton-nucleon scattering 
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e(fci, a) + P{PU V) e'ih) + p'{p2) (1) 

{kl = kl = rn?, Pi = P2 = can be described by the following three 
variables: 



= -q'' = -ih - k2)^ S = 2km, 0, (2) 

where is the azimuthal angle between the scattering plane (ki,k2) and the 
ground level. The Lab system is used with OZ axis along the beam direction 
and plane OZX parallel to the ground level. The explicit expressions of polar- 
ization vectors {C,l and 77) and four-momenta reconstructed in the lab system 



are presented in Appendix A 



The description of the phase space of the radiative process (Fig{T](d,e)) 

e(fci, ^l) + piPi, V) e'ik2) + p'{p2) + 7(fc), (3) 

(fc^ = 0) requires three new kinematic variables: a virtual proton transfer 
momentum squared t = —{ki — ^2 — kY, the inelasticity v = {p2 + k)"^ — M^, 
and the azimuthal angle 0^ between the planes (q, k) and (ki,k2). This set 
of variables defines the four-momenta of all final particles. 

The simulation of radiative events requires an additional definition of the 
lowest bound of the photon energy (or another respective quantity, inelasticity 
Vmin in our case) separating the photon phase space into the region of soft and 
hard photons. Only hard photons need to be simulated while soft photons 
cannot be simulated because of the infrared divergence. The observed cross 
section can be presented in terms of two positively definite parts: 

Cobs = (^rad{Vmin) + BSV (f^min) ■ (4) 



The first term, <Jrad{vmin), describes the cross section with an additional hard 
photon emitted, and the second, crBsv{vmin), contains the contributions of the 
Born cross section, soft-photon emission, and virtual corrections. Here and 
later we define a = da/dQ'^dcp. Note that aobs does not depend on Vmin while 
terms (Tbsv and arad do. 

The strategy for simulation of one event can be defined in a standard way [22]: 

• For the fixed initial energy, Q^, the angle (p, and the missing mass square 
resolution Vmin, the two positively-definite contributions to the observed 
(radiative-corrected) cross section (Jobs, o'radivmin), and (JBSvivmin) are cal- 
culated separately. 

• The corresponding channel of scattering {i.e., BSV or radiative process) is 
simulated for this event in accordance with partial contributions of these 
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two positive parts into the total cross section. More specifically, the channel 
of scattering is simulated in accordance with the Bernoulli trial where the 
probability of "success" {i.e., radiative channel) is calculated as a ratio of 
the radiative part of the cross section to the total cross section. 

• For the radiative event the kinematic variables t, v and (pk are simulated in 
accordance with their calculated distributions. The distributions of v and 
(pk are conditional {e.g. , v is simulated conditionally on t, and (pk is simu- 
lated conditionally on t and v). The explicit expressions for the probability 
densities of these variables are defined by eqs. [T3 

• The four-momenta of all final particles in a required reference frame are 
calculated. 

The initial values of (and (p) can be non-fixed but externally simulated 
according to a probability distribution (for example, the Born cross section). 
If the distribution is simulated over the Born cross section, then the realistic 
observed distribution is calculated as sum of weights computed as ratios 
of the total and Born cross sections for each simulated event. If the observed 
cross section is used for the simulation of Q^, then reweighting is not required. 



3 Explicit expressions for arad{vmin) and aBsv{vmin) 



The analytical expressions for the lowest order RC on which the ELRAD- 
GEN is based, were obtained in ref. [11] (see eqs. (50) and (51)). The 
result for the observed cross section can be formally outlined as aobs = 

/dv 
— [(yR{v) — (To], where C is a kinematic coefficient propor- 
V 

tional to a and the quantity CTji{v) is proportional to the bremsstrahlung 
cross section {(Tr{Q) = (To). This expression does not reproduce the form of eq. 
(jlj), because the term with the integral is not positively definite and the term 
with aR{v) cannot be separated because it is singular for w — >■ 0. Instead, the 
following transformation of this term was used: 



/ —{(^r{v) - o-q) = / —aR{v)0{v - Vr, 

J V J V 



dv f dv 

— ao9{v - Vmin) + J —[crR{v) - Cro]0{Vmin " v). (5) 

The first term in represents the contribution of hard photons, i.e., with in- 
elasticity above Vmin- This term is positively definite and it is used as <Jrad{vmin) 
in (jl]). Its structure and explicit expressions are discussed in Section I3l2l The 
second term admits the analytic integration resulting in correction ^"''^'^{vmin)- 
This term (as well as the third term in the eq. (j5]) discussed in Section 13. 3p 
contributes to the (JBSv{vmin) that represents the part of the observed cross 
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section not contained in the contributions of radiated photons with inelasticity 
above v^,i„ . 



3.1 BSV cross section 



The i^SV-part of observed cross section includes the Born cross section 

(FiglH (a)), loop effects (Fig{T] (b,c)) and the contribution of soft photons. 
The latter is restricted by the inelasticity value v < Vmin'- 

+af{v^,n). (6) 
The Born contribution to the cross section reads: 

^o = ;^E^f-F.(g^). (7) 



The kinematic coefficients 9b are presented in Appendix B The structure 
functions Fi are the squared combinations of the electric and magnetic elastic 
form factors: 



l + Tj 



Ge{Q'')-Gm{Q'') 



"P 



•^4(g') = -m'Gm{q 

with Tp = Q^AM^. 

The factorizing corrections in the first term of ([6]) describe the effects of loops 
and soft-photon emission. The correction S^f comes from the emission of soft 
photons, the 6vr appears as a result of an infrared cancellation of real (Fig{T] 
(d,e)) and virtual (FigJT](c)) photon contribution. The explicit expressions for 
them are: 



Jin f 

TT 



m2 J ^S{S-Q^] 



>VR '■ 



log ^ - 2 - - log - — — + Li2 1 



'k\2 ^m? 2 ^ S-Q^ V S{S-Q^' 

r2- 



IT 



(9) 
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Fig. 2. The region of integration over v- and t- variables for JLab kinematics (Q^ = 3 
GeV^, S = 7.5 GeV^). The hne v = Vmin sphts it into hard (sohd hnes) and soft 
(dashed hnes) real photon regions. 

where Li2 is the Spence function. 

The effect of vacuum polarization by leptons (hadrons) depicted on Fig{T] (b) 
is described by 6^^^ i^vac)- The explicit expression for 5^^^ is defined by eq. (21) 
of ref. [5] while the fit for S!^^^ has been taken from [23] . 

The term S°''^'^{Vmin)co in the R.H.S. of (|6]) contains the correction coming from 
the second term in (ISl). 



6^''{v^^n) = --(log % - l) log (10) 

The last term in (El) is discussed in Section [3l3l 



3.2 Bremsstrahlung cross section 



Since the structure functions depend only on t, and therefore integrals over 
other variables {i.e., v and (pk) can be evaluated analytically or numerically 
with high precision, a reasonable sequence of integration variables is chosen 
such that integration over t is external. This approach allows us to speed up 
the generation of radiative events. The radiative photon phase space for t- and 
f-variables are presented on Fig. [2l It is separated into hard and soft photon 
emission by the line v = Vmin- The cross section of hard-photon bremsstrahlung 



8 



IS 



a™'^(t^™n) = j dtJ2^9f{v„v^a.). (11) 



The quantities 9f-{vi,Vmax) result from the integration over inelasticity v. 
Their arguments correspond to the limits of integration: 



ef{v,,V^a.)=t: J dvW-'e^^iv). (12) 



Here R = + v — t, and the upper sum limits are defined as fcj = (3, 3, 4, 5). 
Accordingly, the quantities OfAv) result from the integration over 0^: 



27r 

ef^{v) = Jd<p,ej'^{v,<Pk). (13) 



The set of quantities 6^ is defined in Appendix B 



Kinematical bounds are defined as 



Jt-Q^)iVi-V4:M^ + t) (t-g')(v^+ V4M2 + t) 
vi = max\ p , p . 



'1,2 



2(M2 + 



(14) 



The probability distributions used for simulation of the photonic variables are 
obtained using ( ITT]) and ( IT^ : 



9 



^^v j=i j=i j=i 
i=l j=l ^^<Afc i=l j=l 



^.c? Contribution of cr]^'^ (vmin) 



The contribution of (j'^'^{vmin) can be presented as an integral over the soft- 
photon region in Fig. |2) 



„addf \ 



AttS^ 



dt / dv 



.i=2 



t2 



(16) 



The hmits of integration over variables t and v read: 



r(t-Q')(v^- V4M2 + t) (t-Q2)(yt + ^4M2 + t) 
fi=max{ p , p }, 

2(M2 + ■ 



The infrared divergences could occur in the limit f ^ — )■ {i.e. at t — )■ Q^) in the 
terms containing . However, one can see that cr 



add( 
R 



is infrared-free. 



Indeed, taking into account Eqs. ( !T3|l . (lB.2p . ( IB.Sp . and (IB.6p . in the limit 



t;^ — 7- we have 



This cancels one degree of R. The second degree of R cancels because the 
integration region is collapsed into a point within this limit. 
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Fig. 3. The structure of the program ELRADGEN 2.0 



4 The structure of the program and input-output data 



4-1 The structure of the program 



The set of files included in the package ELRADGEN 2.0 contains three 
FORTRAN files (eiradgen.f, run.f, test.f), six INCLUDE files (const. inc, 
grid.inc, output. inc, par. inc, pol.inc, test. inc), two data files (rnd.dat, 
test.dat), and one Makefile. No installation is required for this code. 

The eiradgen.f is a source code of the Monte Carlo generator ELRADGEN 
2.0. It contains the set of functions and subroutines for simulation of a single 
event. The loop over simulated events as well as initialization of constants re- 
quires coding in the external program. Two versions of such external programs 
are given in files run.f and test.f. The file run.f is a typical external program 
for simulation of an event with fixed and 0. The file test.f is designed to 
run several tests discussed below. 



The structure of the code is illustrated in Fig. El 
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• program main is a sample of an external program that invokes ELRAD- 
GEN (in our case it is the program included in run.f and test.f); 

• elradjnit defines all constants (such as beam energy and polarization 
degrees) which are necessary for generation; 

• gridJnit prepares the grids for generation of photonic kinematic variables; 

• eiradgen is a main subroutine governing the simulation of an event; 

• urand is a generator of uniformly distributed random numbers; 

• fsib calculates the Born cross section; 

• fsigt invokes one of the subroutines fsigtan or fsigtnu to calculate the 
cross section dcr/dt; 

• fsigtan calculates the analytical cross section da/dt for unpolarized scat- 
tering; 

• fsigtnu calculates the cross section da/dt with numerical integration over 
variable v for polarized scattering; 

• fsigt V calculates the analytical cross sections da/{dtdv) and da/{dtdvd(l)k); 

• ffpro is a model for elastic form factors. 

The six INCLUDE files are: 



const. inc includes all necessary constants, e.g. , the fine electromagnetic 
constant, the proton and lepton masses; 

grid. inc includes nets of bins for simulation of the three photonic variables; 
output. inc contains variables governing the form of output as discussed 
below 

pol.inc includes quantities which describe the polarization state (defined 
in 



Appendix A.l ) 



• test. inc includes variables and nets of bins required for test run; 

• par. inc includes variables required for calculation of cr'^'^ivmin)- 

The file rnd.dat includes an initial integer for the flat generator urand, and 
test.dat is an example of output data for a test run (when test.f is used as an 



external program); the results of different test are presented in Appendix C 



The commands "make" or "make test" need to be run for creating the exe- 
cutable file for the simulation or for the test runs, respectively. 



4.2 Input- output data 



Input data in ELRADGEN 2.0 are set up in program main of run.f or test.f. 
Majority of them are transferred to the main program through parameters in 
the subroutine eiradgen. They are: 

• ebeam is an energy of electron beam; 

• q2 is a virtual photon momentum squared Q^] 
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• phi is an azimutlial angle between tlie scattering plane and the ground 
level; 

• vvmin is a missing mass square resolution Vmin for separation of radiatively 
corrected cross section into radiative and BSV parts; 

• vcut is a cut-off quantity Vcut that allows to exclude the simulation of hard 
photons above Vcut- 

The last variable provides the opportunity to exclude simulation of events 
with inelasticity above a predetermined level. This could be convenient when 
simulation is performed for experimental design, when hard real photon are 
removed from experimental data by putting a cut on the missing mass of the 
undetectable particle. 

The quantities describing the polarization characteristics of beam and target 



(defined in Appendix A.l ) are transferred to the code through the common 
block pol containing four variables: i,ii) plrun and pnrun, the polarization 
degrees of the lepton beam and target Pn, iii) thetapn, the angle 9r^ between 
3- vectors of the target polarization r] and initial lepton momentum ki, and 
iv) phipn, the angle 0^ between OZX and (ki, rj) planes. 

One additional variable itest governs the form of the output. If itest 7^ 0, all 
output information is printed to the file test.dat. If itest = 0, the output 
data are collected in two common blocks of the file output. inc: 

common / variables / tgen,vgen,phigen,weight ,ich 



and 



common / vectors / vprad,phrad 

Here tgen, vgen, and phigen are the generated photonic variables t, v, and 
(pk, respectively, weight is a ratio of the observable cross section to the Born 
one, variable ich shows whether the scattering channel is radiative {ich = 1) 
or BSV {ich = 0). The quantities vprad = P2 — Pi and phrad := k are four- 
momenta of virtual and real photons defined in the Lab system. 

For BSV events, vgen = 0, tgen = Q^, (pk = 0, phrad = and vprad = k\ — k2. 



5 Numerical tests and comparison with experimental data 



Below we describe three types of numerical experiments allowing: i) to cross- 
check some key distributions and parameter estimates in ELRADGEN, ii) to 
investigate issues related to a possible dependence of simulated cross sections 
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Fig. 4. Histogram (points) and corresponding probability densities (solid lines) for 
variables describing the exclusive real hard photon production in polarized elec- 
tron proton scattering at JLab kinematic conditions {Ebeam = 4 GeV, = 3 
GeV^, ) for transverse polarized proton {6rj = 90^) with (p = cj)^, PlPn = —1 and 

Vrmn = l^'^ GeV^. 

on Vmini and iii) to perform a comparison with data collected in the BLAST 
experiment. 

5.1 Tests implemented in ELRADGEN 

There are five tests implemented in the program. The first three deal with 
checking to which extent the simulated distributions on photonic variables t, 
f , and (pk correspond to analytical probability distributions given in eq. (fTSll . 
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Fig. m presents the t-, v-, and 0fc-distributions calculated numerically and gen- 
erated by ELRADGEN under JLab kinematic conditions for a transverse po- 
larized target. The theoretical and simulated distributions of this case (as well 
as for unpolarized and longitudinally polarized targets) are almost identical. 

The sharp peaks in the t-distribution coming from the collinear singularities, 
i.e., from the kinematical region where the real photon is emitted along either 
the initial or the final lepton. After integration over the inelasticity v, these 
two singularities are situated near t = and are only slightly different. 

The peaks on the plots of the f-distributions correspond to the collinear singu- 
larities as well. Since the variable t is external, the f-distribution is conditional 
on t, and therefore only one peak corresponding to either the initial or the 
final electron appears for each f-distribution. 

Finally, the ^^-distributions show that most of the photons are emitted in the 
scattering plane. 

For generation of the three types of distributions, one has to set itest = 1, 
itest = 2, or itest = 3, in the file test.f, and then type "make test" and 



"./test.exe". In Appendix C, the test outputs for t, v, and (pk generation with 
P = 1, i?^^^=4 GeV, 6rf = 48°, 20 bins for the histogramming and 10^ radiative 
events are presented. 

After generating all photonic variables for one radiative event, ELRADGEN 
reconstructs the four-momenta of the final particles. To make sure that the 
vectors are constructed properly, the next test corresponding to itest := 4 
is implemented. This test allows to perform the numerical comparison of the 
generated variables t, v, and (j)k with the value of these variables reconstructed 
from four-momenta of the particles. This test also reconstructs the mass of 
the real photon that has to be equal to zero. 

The test with itest := 5 provides us with the comparison of the unpolarized 
cross section integrated over v analytically and numerically. 



5-2 Vmin- dependence and comparison with MASCARAD 



The Monte Carlo generator ELRADGEN 2.0 was developed on the basis of 
the FORTRAN code MASCARAD, therefore the agreement of outputs of 
both programs with the same input parameters has to be demonstrated as a 
primary test. Here we restrict our crosscheck to the JLab kinematic conditions 
without cuts on inelasticity v and focus on the comparison of the ratio of the 
radiatively corrected cross section to the Born one. Define components of the 
cross sections as: 
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Table 1 

The ^^mm"dependence of the ratios of radiative, BSV, and observable contributions 
to the unpolarized {PlPn = 0) electron-proton cross section to the Born cross 
section for JLab kinematic conditions {Ebeam = 4 GeV and = 3 GeV^): a (b) 
presents the results of analytical (numerical) integration over v in ELRADGEN, 
while c shows the results of the calculation using MASCARAD |11] . 
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Table 2 

^'mm-dependence of the radiative, BSV, and observable contribution to electron- 
proton scattering with longitudinally polarized target {Orj = 0) for a different spin 
orientation in the Born units and comparison with MASCARAD [llj for JLab 
kinematic conditions (-Ebeam = 4 GeV and = 3 GeV2). 

a^'^^i^L, VL,t) = < + PLPNCr'ai^L, VL,t), (19) 

where a = 0, BSV, rad, obs 

First, we consider the unpolarized scattering for which an option with an- 
alytical integration over v is available. Table 1 presents the results of the 
analytical and numerical integration for the BVS- and radiative contributions 
{i.e., (iBsvivmin) and (Tradivmin)) to the observed cross section, as well as the 
results obtained using MASCARAD. Each of these contributions changes es- 
sentially by decreasing Vmin from 1 to IQ-^ GeV2, while the observable cross 
sections barely change for both analytical and numerical integration over v. 



16 



(jiG V 


T 1 T 


T / T 


weight = cr^^Ja^ 


ELRADGEN 


ELRADGEN 


ELRADGEN 


MASCARAD 


TD ID 


i 


-i 


i 


-i 


i 


-i 


i 


-i 


i 


0.9457 


1.447 


0.9730 


0.9746 


1.919 


2.422 


1.917 


2.420 


in-i 


1.117 


1.620 


0.8018 


0.8018 


1.918 


2.422 


10-2 


1.273 


1.776 


0.6386 


0.6386 


1.912 


2.415 


10-3 


1.432 


1.935 


0.4754 


0.4754 


1.908 


2.411 


10-4 


1.596 


2.099 


0.3122 


0.3122 


1.908 


2.411 



Table 3 

■Wmm-dependence of the radiative, BSV and observable contribution to electron- 
proton scattering with transversely polarized target (0^ = '7r/2, (j) = 0^) for different 
spin orientation in the Born units and comparison with MASCARAD [TT] for JLab 
kinematic conditions (-Efceam = 4 GeV and = 3 GeV^). 

The similar behavior of the radiative and BSV parts takes place for the polar- 
ized case. This is illustrated in Tables 2 and 3. The observable cross sections 
change by no more than 1%. 

5.3 Results from the BLAST data 

Analyzing the A-excitation region in ep-scattering, it is necessary to extract 
the contribution of real hard-photon emission that accompanies the elastic 
ep-scattering (so-called elastic radiative tail) and cannot be removed from the 
data by any experimental cuts. The main radiative photons are emitted by the 
electron leg (see FigJT](d,e)), because their contributions include the logarithm 
of the electron mass. These radiative events are spin-dependent, and therefore 
affect not only the cross section, but other extracted quantities in the A- 
region as well, e.g. asymmetries, spin-correlation parameters, spin-structure 
functions, etc. 

The BLAST experiment was designed to study spin-dependent electron scat- 
tering off protons and deuterons with small systematic uncertainties [25]. The 
experiment used a longitudinally polarized, an intense electron beam and 
isotropically pure highly-polarized internal targets of hydrogen and deuterium 
from an atomic beam source. For extraction of the elastic radiative tail con- 
tribution, the new version 2.0 of Monte Carlo generator ELRADGEN has 
been applied. This generator was incorporated into the BLAST Monte Carlo 
event generator [26], where longitudinally polarized electrons at an energy of 
850 MeV and at a polarization factor of 65%, were scattered off a highly- 
polarized hydrogen internal gas target (Pat ~ 80%), with the average target 
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Fig. 5. Normalized yields as a function of the invariant mass, Vl^[GeV] over 
0.08 < Q2 < 0.38 GeV^ The dots show the BLAST ABS hydrogen data corrected 
for the background contributions, and the solid line represents the Monte Carlo 
simulations with radiative effects (ELRADGEN 2.0). 




Fig. 6. The effect of the radiative contributions to the asymmetry in the A-excitation 
region. The left (left) and right (right) asymmetries are shown with (dots) and 
without (squares) radiative corrections (RC), for 0.08 < < 0.38 GeV^. Monte 
Carlo simulations using the MAID 2003 model [27] (straight line) and ELRADGEN 
(dotted line) are shown for comparison. 

spin direction oriented at 48.84° to the left of the beam direction. 
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In order to estimate the contribution of the elastic radiative tail to the A- 
excitation region, the results of the Monte Carlo simulations were normalized 
to the data elastic peak, as shown in Fig. [51 In this figure, the normalized yields 
correspond to the outgoing electrons detected in each sector of the BLAST 
detector (inclusive scattering). The radiative tail obtained from the above 
normalization is subtracted from the measured yields (radiative corrections), 
and then the left and right asymmetries are extracted. For comparison, in 
Fig. |6] we show the left and right asymmetries with and without the radiative 
corrections. The asymmetry from ELRADGEN alone is also shown (dotted 
line) in order to see the radiative tail effect to the overall asymmetries and its 
spin dependence. 



6 Conclusion 

In this paper, we presented a new version of Monte Carlo generator ELRAD- 
GEN for simulation of real-photon events within elastic lepton nucleon scat- 
tering for longitudinally polarized lepton and arbitrary polarized target. Fol- 
lowing the absolute necessity of both accuracy and quickness for our program, 
we have developed the fast and highly precise code using analytical integra- 
tion wherever it was possible. The developed program has a broad spectrum 
of applications in data analysis of various experimental designs on polarized 
ep-scattering, including the measurements of the generalized parton distribu- 
tions, the generalized polarizabilities, and the evaluation of spin asymmetries 
in elastic scattering. Also, it can be used as a generator of the "Born" process 
in DVCS measurements and of the radiative tail from the elastic peak in DIS. 
The most significant application of the generator is in the experiments with 
the complex detector geometry. 

The set of numerical tests of the presented version of this code proved its 
high quality. First, a good agreement with FORTRAN code MASCARAD [11] 
was found. Second, no dependence on the missing mass square resolution was 
found. Third, the distributions of the generated radiative events are found 
to be in accordance with the corresponding probability densities. Fourth, a 
good agreement with the radiative tail from the elastic peak measured in the 
BLAST experiment was demonstrated. 

Several additional steps allowing to make the simulation even faster are 
planned. They include implementation of analytical integration over variable 
V of the hard-photon emission contribution with the longitudinally polarized 
target and utilizing the look-up table option for faster simulation of radiative 
events with transverse component of the target polarization vector. 
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Fig. 7. 3-vectors decomposition in the LAB system 

The spectrum of apphcations of the presented code could be extended after a 
certain substantive upgrade in several directions including the development of 
this generator for transferred polarization from lepton beam to recoil proton 
[28] , and for involving this generator into the measurement of the electromag- 
netic form-factors of the proton in elastic scattering with unpolarized [29] and 
polarized targets [30)131] . Inclusion of electroweak effects will provide the gen- 
eralization for the investigation of electroweak corrections in experiments on 
axial form factors of the nucleon [32] and parity violation elastic scattering 
[33] ■ This generator can be included in data analysis of experiments with the 
measurement of unpolarized and spin-flip generalized polarizabilities in virtual 
Compton scattering [?lfi] . 
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Appendix A Four-vectors 

In this section we present the explicit expression for four-vectors decomposi- 
tion in LAB system depicted in Fig. [7] 
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Appendix A.l Polarization vector definitions 



As it was mentioned above we assume that the electron beam has a longitu- 
dinal polarization. Therefore its polarization vector has a form 



iL = ^{-h-2mpA. (A.l) 



The target polarization vector in the Lab. system can be decomposed into 
longitudinal 



VL = ^ (2Mh - ) (A.2) 



V ' M 

and transverse r/y components as it is depicted in Fig. [7] 

ri = cos{9^)r]L + sm{9r^)r]T, (A. 3) 

where 6r^ is the angle between 3-vectors ki and 77. Transverse tjt component 
can be presented as: 

rjT = cos(0 - (f)^)T]t + sin(0 - (f)r,)v±, (A.4) 
where 0.^ is the angle between (ki, rj) and OZX planes, and 

(4m2M2 + - SX)ki + A.fca - (SQ^ + 2m^S^)pi 



\ = (a , 

|k2||ki| sin 6' 



ry^= (0, ^^^^^^^ A = (O,sin0,-cos0,O). (A.5) 



Here 



X = 2k2Pi, S,^ = S-X, A = SXQ^ - M^Q^ - m^Ag, 

\g = Sl + 4Q^M^. (A.6) 

It should be noted, that for the BSV process the variable X is fixed by 
and S: X = S — Q^, while for the radiative one, the variable X (as well as 
and Xq) depends on inelasticity: 

X = S - - V, = + V, Xg = {Q^ + v) + AQ^M\ (A.7) 
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As it follows from f lA.9P the normal to scattering plane component of t] satisfies 
the equations: 



kiri± = k2'n± = PiV± = 0, kr]± = -p2'n± = sincpk- 



(A.8) 



Appendix A. 2 Four-momenta reconstruction 



After generation of photonic variable t, v and (pk the four-momenta of final 
proton p2 = {p^2\p2 ■,P2'\p2^)^ lepton ^2 = (fc2°\ ^2^^ ^2^^ M^'') ^-iid real pho- 
ton k = (A;(°\ k'^^\ k^^\ A;^^^) in the LAB system read: 



(2) 
P2 



(3) 
P2 



rvn 



'A3( Ai cos cos 0fc — ^JKiS sin sin + A2V^ cos ( 



^^{\fKiS cos (f) sin + Ai sin cos 0^) + A2-\/A4sin0 

t + 2M^ 



A1A2 - 4M^ VA3A4 cos 0fc (0) 

2AgM^ ' P2 - 2M ' 

/Alsin0 (3) _ - Ai (0) _ S -Q^ - V 

7; } '^2 — r^A ) '^2 



/AlcOS0 (2) _ _ 

S ' ^2 " S '"^ ~ 2MS ' 



2M 

^^{^JX^S sin (j) sin (f)k — Xi cos cos 0^) + (A^ — A2)a/A4Cos0 

XqS 

— \^{\f\S cos (f) sin (f)k + Ai sin cos 0^) + {Xq — A2)-\/A4sin( 

AgS 



Ai(Ag - A2) + 4M^v/A;A;cos0fc ^ + - t 
2AqM5 ' 2M 



(A.9) 



where 



Ai = S{Q^ + v) + 2M^Q\ X2 = t{Q^ + v) + 2M\Q^ + t), 

A3 = tv{Q^ -t + v)- M\Q^ - t)\ A4 = Q^S{vmax ' v) . (A.IO) 



Appendix B Explicit expressions for the kinematic quantities 6 



The kinematic coefficients Of appear as a convolution of the leptonic tensor 
with corresponding hadronic structures 
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IIV 

Wi 



-iPive 



M 



111/ 

W2 



W 



P^^Pu 



(B.l) 



and read 











oi = 









{S{S - Q'') - M^Q^' 



2M2 



PlPn- 



^^^^^^ 
mQ"^ qr] 
M3 



(B.2) 



Here Pl and Pjv define the degree of the lepton and target polarization re- 
spectively and the explicit expressions for polarized vectors of the scattering 



particles ^ and r] can be found in Appendix A.l 



The quantities Oi{vi,V2) appear as a convolution of the leptonic tensor that is 
responsible for the real photon emission: 





^kki 


k2a^ 

kk2 J 


7m- 


1 ilk-fa 

2kki 


2kk2 




^kia 
[kki 


k2a^ 

kk2 J 




lak^y 

2kki 


li/kja 

2kk2 



(B.3) 



with hadronic structures presented in eq. (lB.ip . 
As a result 
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Appendix B.l Quantities 6 ij{v) and 9ij{v, 



Here, we combine explicit expressions for 0^{v) and 9^{v, (pk) quantities cal- 
culated for polarized scattering j5|19|34] and present the explicit expressions 
Of-{yi,V2) calculated for unpolarized scattering only. 

Both types of quantities 6^j{v) and 9fj{v,(t>k) for i = 1,2,3 take the similar 
form 



e^, = 2{SX - M^Q^)Fir/M^, 

eg = {2m'SpF2. + S,S,F^+ + 2{S, - 2M^t)Fir - tSIF,)/2M\ 
eg = {4M^F + (2MV - S,)tF, - SpF,+)/2M\ 

eg = PiPN^im hi - iv)FiR, 

eg = -P^P^^(2r/g {ThiF, - 2F|^) + Q'vIC{fI - Fi, - 2F^) 

+iivrFiR - Am'him - ^2+)) 
Og = PLPN^^ivlC r(2Fi + Fl - F^) - 2 hi tF] - Am'F^'^ - 6Fp^ 

+Q\Fll-Fl^_)l 

eg = PlPn^ ( 2^1" + Fll - Ft ) , (B.5) 

where t = {t — Q^)/R and the four- vector }C = ki + k2. 
For i = 4 we have 



ogi = mOAi/M, eg^ = iriqe,2 - e2,)/M, 

og = {m043 - ei2)/M, eg = {mh^ - 01^)Im, 

Og = -e2jM, (B.6) 
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where 



4-77? 

e,i = PlPnj^ { 2ipQ^-S, ) Fir, 



m 



9,3 = PlPnj^HQ' - tS, ) Fl -(Ql-rSp) fI^ + 2 k^i tF^ 

+6F/^-2Fir5,), 
eu = -PlPn^, ( 2F| - F|_ + F|+ ) . 
The quantities 92j are calculated as: 



(B.7) 



92j - 9Aj{Faii F^ii, F^ii F^H). 
The upper indices in Fan appears in the following way: 



(B.8) 



2F|+ = (2Fi+ + rF2_)sj + Fg+r^, 

2^2^ = (2Fi+ + tF2-)s^ + F2+r^ + 

2F|_ = (2Fd + F2+)rs^ + F^.r^, 
2F| = Fi+s^ + F,r^, 



4 sin 



^ A, 



A., 



'2+, 



2F2 = F^+Sr, + Farr,+ 



4 sin 



4F|^ = (2Fi+ + rF2_)(r^sj + s^r^) + F2+(r^r5 + r's^s^) 



+4(2F + FrfT=')s^s^ + 



8 sin (t)kdr. 



R \\ 



2+1 



8 sin 0fe<^?7 



^ A, 



\ 2-) 
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" sin (pi^dfj 



+ 



The quantities 



R 



(B.9) 



^{(,v} = ^{(M + ^{^M'^i^M = '^i^UM ~ ^{i,ri}) + 2c{c,77} (B.IO) 
are combinations of coefficients of polarization vectors ^ and 1] expansion over 



basis (see [Appendix A.T]) 



^ = 2(a^fci + b^k2 + c^p), 



(B.ll) 



We note that the scalar products from (lB.5fB.6fB.7p are also calculated in 
terms of the polarization vector coefficients: 



-^7] = 2m {a^a^j + 6^6^) + 2M c^Cr, + Q^(ag6,, + b^a^) 

+S{a^Cr^ + C^ttrj) + X{b^Crj + c^b^). 



(B.12) 



With the exception of the contribution proportional to drj all dependencies of 
9ij on the photonic variable (pk are included in the quantities F, however in 
both cases we have: 



Fir = m^F2+ - Q^F^ 



(B.13) 



So for %(f , (pk) the quantities F read 



Fd{v,(pk) 



F{v, 



Z\Z2 



Fi+{v,(pk) = F{v,(pk) 



F2±{v, (pk) = F{v, (pk) {\±\] , F{v, (pk) 

\Z2 4 



2txJ\ 



1 1 

Zl Z2 



(B.14) 



Here 
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2kki 1 I — 

K \q 

2kk'i 1 I — 
Z2^^ = —{Q''Sj, + t{XS, - 2M2g2) - 2MJx,cos(j)k): (B.15) 

K Aq 



and 



Sp^ S + X ^2S -Q"^ -V. (B.16) 
The following equalities define the functions F for Ofj{v): 

F2±{v) = B,{t)C:^"\t) t B,{T)Cf'\T), (B.17) 
where 

Si,2(t) = -\ ( XqT ± 5p(5,T + 2Q^) ) , 

C,(r) = {St + g2)2 + 4m2(g2 + - t'M^), 

C2(t) = (Xr - g2)2 + W(g2 + tS-, - r'M^). (B.18) 

We note that F^ has a 0/0- like uncertainty for r = (inside the integration 
region). It leads to difficulties in numerical integration, so another form is used 
also 



Appendix B.2 Quantities 9^{vi, V2) 

As it was mentioned above for the unpolarized scattering, the integration over 
V is performed analytically resulting in: 

(^1, V2) = -4If + Atl^^ - 2{Q^ + 
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0^{vu V2) = ^("^M'If - t{ll + 2J°) + t{2S - t)I^l + < 

-4(25 - t)/2i^ + 4(S^ - t{M^ + S))I^l 
+t{Q^ + AS- 3t)Ia^ 

+[tiQH - {2S - tf) + 2M\t' + Q')]IJ^), (B.20) 

where 



1° 
-'21 



: / dvFiv) = log ^ I = 

J \ n2 _L _L ./('/n2 I „ \2 I 4A/f2n 



V2 



«2 



1 ? 1 92 

= -m' j dv{F,4v) - F,.{v)) = -iQ' - t)-jA?, 



Vl 
V2 



7^^ 



' t3 2-^2(g2_^) Ige ^3 I' 



.1 Idv^, ^ S S-t 1 /A} A^' 



1)2 



^AL.-iAL. . (B.21) 



Vl 



Here 



A2 = |A(^2)|AM - |A(^i)|A(^i), = \D,{V2)\ - mv^)\, 
A° = D,{v2)/Dl,{v2) - Di{v{)/Dl^^{v,), (B.22) 
r 2m2t(g2 - t + 2^2) + Q^Di{v2) + A(t^2)\/4m2t + 
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log 



'2m^Vi{Q^ - t + 2v2) + ViD2{v2) + 1^5(^2) V4m2 + t 
_2m^Vi{Q^ - t + 2vi) + ViD2{vi) + Dsl^^i) V4m2 + t_ 



and 



D,{v) = {t-Q^){S-Q^) + Q\, 

D2{v) = S{Q^ -t)+tv, 

Dsiv) = vt{v -t + Q^)- M^Q^ - t)\ 

D,{v) = ^Dl{v)+AmWsiv). (B.23) 



Appendix C Test output 



Here, we present the results of the test as test.dat output file corresponding 
to: 

1) itest := 1 - the generation of p{t) distribution and comparison with the 
analytical cross section corresponding to the first formula in (|T5|) (here and 
below invariants v,t are in GeV^) 

itest=l 

t generation 

rgen is generated probability 

rcalc is calculated probability 

Ebeaiii=.850 GeV 

Q**2=.200 GeV**2 

vinin=0. lOE-01 GeV**2 

vcut=0.00 GeV**2 

PL*PN=-1.00 lepton polarization times nucleon polarization 
thetapn=48.0 degrees angle between target polarization 

vector and beam momentum in Lab system 
phipn=0.00 degrees athimutal angle in Lab system 
number of bins 20 

number of radiative events 100000000 
initial random number 12 

bin t rgen rcalc rgen/rcalc 

GeV**2 GeV** (-2) 

1 0.4237E-01 1.624 1.786 0.9093 

2 0.9197E-01 1.619 1.589 1.019 

3 0.1453 2.325 2.274 1.023 

4 0.1969 12.14 35.10 0.3460 

5 0.2351 1.518 1.666 0.9108 



29 



6 





.2893 


0, 


.3042 




0. 


.3153 




0.9647 


7 





.3411 


0, 


.1113 




0. 


,1134 




0.9816 


8 





.3923 


0, 


.5102E- 


-01 


0. 


.5163E- 


-01 


0.9883 


9 





.4433 


0, 


.2666E- 


-01 


0. 


.2688E- 


-01 


0.9917 


10 


0, 


.4943 


0. 


, 1518E- 


-01 


0. 


, 1526E- 


-01 


0.9945 


11 


0, 


.5449 


0. 


.9216E- 


-02 


0. 


,9259E- 


-02 


0.9954 


12 


0, 


.5960 


0. 


,5872E- 


-02 


0. 


,5865E- 


-02 


1.001 


13 


0, 


.6466 


0. 


, 3893E- 


-02 


0. 


,3872E- 


-02 


1.006 


14 


0, 


.6973 


0. 


,2577E- 


-02 


0. 


,2632E- 


-02 


0.9791 


15 


0, 


.7479 


0. 


. 1822E- 


-02 


0. 


, 1833E- 


-02 


0.9937 


16 


0, 


.7983 


0. 


. 1278E- 


-02 


0. 


. 1304E- 


-02 


0.9801 


17 


0, 


.8490 


0. 


.9465E- 


-03 


0. 


.9396E- 


-03 


1.007 


18 


0, 


.9003 


0. 


.6965E- 


-03 


0. 


.6806E- 


-03 


1.023 


19 


0, 


.9509 


0. 


.4965E- 


-03 


0. 


.4945E- 


-03 


1.004 


20 


0, 


.9973 


0. 


.2696E- 


-03 


0. 


.3558E- 


-03 


0.7577 



2) itest := 2 - the generation of p{v) distribution and comparison with the 
analytical cross section corresponding to the second formula in f|T5l) 

itest=2 

V generation 

rgen is generated probability 

rcalc is calculated probability 

Ebeam= . 850 GeV 

Q**2=.200 GeV**2 

vmin=0. lOE-01 GeV**2 

vcut=0.00 GeV**2 

PL*PN=-1.00 lepton polarization times nucleon polarization 
thetapn=48.0 degrees angle between target polarization 

vector and beam momentum in Lab system 
phipn=0.00 degrees athimutal angle in Lab system 
number of bins 20 

number of radiative events 100000000 
initial random number 12 

t= 0.5218 GeV**2 
bin V rgen rcalc rgen/rcalc 





GeV**2 




GeV**(- 


-2) 






1 


0.6260 


0, 


.1660 


0. 


,1657 


1.002 


2 


0.6599 


0, 


.1931 


0. 


,1933 


0.9990 


3 


0.6937 


0, 


.2265 


0. 


,2266 


0.9998 


4 


0.7275 


0, 


.2676 


0. 


,2677 


0.9996 


5 


0.7613 


0, 


.3206 


0. 


,3205 


1.000 


6 


0.7952 


0, 


.3904 


0. 


,3913 


0.9978 


7 


0.8291 


0, 


.4927 


0. 


,4921 


1.001 


8 


0.8630 


0, 


,6478 


0, 


,6486 


0.9987 



30 



9 


0.8971 


. 9284 


0.9278 


1.001 


10 


0.9316 


1.579 


1.580 


0.9993 


11 


0.9698 


6.013 


6.043 


0.9951 


12 


0.9886 


13.45 


16.44 


0.8179 


13 


1.029 


1.759 


1.760 


0.9992 


14 


1.064 


0.9512 


0.9502 


1.001 


15 


1.098 


. 6343 


0.6361 


0.9973 


16 


1. 132 


0.4694 


0.4686 


1.002 


17 


1.166 


0.3649 


0.3645 


1.001 


18 


1.199 


0.2937 


0.2942 


0.9985 


19 


1.233 


. 2435 


. 2436 


0.9999 


20 


1.267 


0.2071 


0.2053 


1.008 



3) itest := 3 - the generation of p{(pk) distribution and comparison with the 
analytical cross section corresponding to the third formula in fllSp 

itest=3 

phik generation 

rgen is generated probability 

rcalc is calculated probability 

Ebeam= . 850 GeV 

Q**2=.200 GeV**2 

vmin=0. lOE-01 GeV**2 

vcut=0.00 GeV**2 

PL*PN=-1.00 lepton polarization times nucleon polarization 
thetapn=48.0 degrees angle between target polarization 

vector and beam momentum in Lab system 
phipn=0.00 degrees athimutal angle in Lab system 
number of bins 20 

number of radiative events 100000000 
initial random number 12 

t= 0.5218 GeV**2 

v= 0.7962 GeV**2 
bin phik rgen rcalc rgen/rcalc 





rad 




rad** (- 


-1) 








1 


-2.984 


0, 


. 1890E-01 


0. 


1884E- 


-01 


1.003 


2 


-2.668 


0, 


.2043E-01 


0. 


2036E- 


-01 


1.004 


3 


-2.351 


0, 


.2388E-01 


0. 


2381E- 


-01 


1.003 


4 


-2.035 


0, 


.3031E-01 


0. 


3027E- 


-01 


1.001 


5 


-1.718 


0, 


.4189E-01 


0. 


4197E- 


-01 


0.9980 


6 


-1.402 


0, 


.6354E-01 


0. 


6381E- 


-01 


0.9959 


7 


-1.085 


0, 


.1062 


0. 


1069 




0.9927 


8 


-0.7682 


0, 


.1961 


0. 


1988 




0.9864 


9 


-0.4529 


0, 


.3922 


0. 


4013 




0.9773 


10 


-0.1470 


0, 


.6982 


0. 


7223 




0.9667 
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11 


0.1470 


0. 


.6981 




0. 


.7223 




0.9665 


12 


0.4530 


0. 


.3924 




0. 


.4012 




0.9779 


13 


0.7682 


0. 


.1960 




0. 


.1988 




. 9862 


14 


1.085 


0. 


.1061 




0, 


.1069 




0.9919 


15 


1.402 


0. 


.6351E- 


-01 


0, 


. 6380E- 


-01 


0.9955 


16 


1.718 


0, 


.4192E- 


-01 


0. 


.4197E- 


-01 


0.9988 


17 


2.035 


0, 


.3036E- 


-01 


0. 


.3026E- 


-01 


1.003 


18 


2.351 


0, 


. 2389E- 


-01 


0, 


.2381E- 


-01 


1.003 


19 


2.667 


0. 


. 2036E- 


-01 


0, 


. 2036E- 


-01 


1.000 


20 


2.983 


0. 


. 1888E- 


-01 


0, 


. 1884E- 


-01 


1.002 



4) itest := 4 - the cross-check of the accuracy of the vector reconstruction for 
5 random radiative events 



itest=4 

variable reconstruction 
Ebeam=.850 GeV 
Q**2=.200 GeV**2 
vmin=0.10E-01 GeV**2 



vcut=0 . 00 
PL*PN=-1.00 
thetapn=48 . 



phipn=0 . 00 
number of bins 
number of radiative events 
initial random number 12 



GeV**2 

beam polarization times target polarization 
degrees angle between target polarization 
vector and beam momentum in Lab system 
degrees athimutal angle in Lab system 
20 

5 



event= 1 
test V reconstruction 
v=0.206247E-01 GeV**2 
v=0.206247E-01 GeV**2 



reconstructed v from 4-vector 
generated v 



test t reconstruction 
t=0 . 202340 GeV**2 
t=0 . 202340 GeV**2 



reconstructed t from 4-vecto 
generated t 



m2gamma= 0.157778E-10 GeV**2 real photon mass square 



event= 2 

test V reconstruction 
v=0. 579509 GeV**2 reconstructed v from 4-vector 
v=0. 579509 GeV**2 generated v 
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test t reconstruction 
t=0. 119041 GeV**2 
t=0. 119041 GeV**2 



reconstructed t from 4-vecto 
generated t 



m2gaiiima=-0 . 855924E-08 GeV**2 real photon mass square 



event= 3 

test V reconstruction 

v=0 . 328652E-01 GeV**2 reconstructed v from 4-vector 

v=0.328652E-01 GeV**2 generated v 

test t reconstruction 

t=0. 203009 GeV**2 reconstructed t from 4-vecto 

t=0. 203009 GeV**2 generated t 

m2gamma= 0.249049E-10 GeV**2 real photon mass square 



event= 4 

test V reconstruction 
v=0. 579878 GeV**2 reconstructed v from 4-vector 
v=0. 579878 GeV**2 generated v 

test t reconstruction 

t=0. 117210 GeV**2 reconstructed t from 4-vecto 

t=0. 117210 GeV**2 generated t 

m2gaimna=-0 . 340940E-08 GeV**2 real photon mass square 



event = 5 
test v reconstruction 
v=0.412302E-01 GeV**2 
v=0.412302E-01 GeV**2 



reconstructed v from 4-vector 
generated v 



test t reconstruction 
t=0. 193803 GeV**2 
t=0. 193803 GeV**2 



reconstructed t from 4-vecto 
generated t 



m2gamma= 0.124488E-10 GeV**2 real photon mass square 



4) itest := 5 - the cross-check of the accuracy of the analytical and numerical 
integration over t for unpolarized scattering 
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itest=5 

comparison of the analytical and numerical integration over t 

rnum is the numerically integrated cross section 

ran is the analytically integrated cross section 

Ebeain=4 . 00 GeV 

Q**2=3.00 GeV**2 

vmin=0. lOE-01 GeV**2 

vcut=.200 GeV**2 

PL*PN= 0.00 beam polarization times target polcirization 



lin 




t 


rnum 




ran 




mum/: 




GeV**2 


nbarn*GeV** (■ 


-4)*rad**( 


-1) 




1 


2. 


.340 


0, 


.8162E-06 


0. 


.8179E- 


-06 


0.9979 


2 


2, 


.382 


0, 


.2651E-05 


0, 


.2656E- 


■05 


0.9980 


3 


2, 


.424 


0. 


.4847E-05 


0, 


.4857E- 


■05 


. 9980 


4 


2, 


.466 


0. 


. 7560E-05 


0, 


.7575E- 


■05 


0.9981 


5 


2. 


.509 


0, 


. 1102E-04 


0. 


. 1104E- 


-04 


0.9981 


6 


2. 


.551 


0, 


. 1561E-04 


0. 


. 1564E- 


-04 


0.9982 


7 


2. 


.593 


0, 


.2194E-04 


0. 


.2198E- 


■04 


0.9983 


8 


2, 


.635 


0. 


.3111E-04 


0, 


.3116E- 


■04 


. 9984 


9 


2, 


.677 


0. 


.4521E-04 


0, 


. 4528E- 


■04 


. 9985 


10 


2. 


.720 


0. 


. 6865E-04 


0. 


. 6874E- 


■04 


0.9986 


11 


2. 


.762 


0, 


. 1122E-03 


0. 


. 1124E- 


-03 


0.9988 


12 


2. 


.804 


0, 


.2096E-03 


0. 


.2098E- 


-03 


0.9990 


13 


2. 


.846 


0, 


. 5433E-03 


0. 


. 5436E- 


-03 


0.9994 


14 


2, 


.889 


0. 


. 5977E-02 


0, 


.5971E- 


■02 


1.001 


15 


2, 


.931 


0. 


.1018E-01 


0, 


. 1017E- 


■01 


1.001 


16 


2. 


.973 


0, 


.2645E-01 


0. 


.2641E- 


■01 


1.002 


17 


3. 


.015 


0, 


.4498E-01 


0. 


.4490E- 


-01 


1.002 


18 


3. 


.057 


0. 


.1006E-01 


0. 


. 1004E- 


-01 


1.001 


19 


3, 


.100 


0. 


.3106E-03 


0, 


.3109E- 


■03 


0.9991 


20 


3, 


.142 


0. 


. 2593E-04 


0, 


. 2598E- 


■04 


. 9982 
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